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1. Introduction 

The main target of lattice studies of tiie Higgs-Yukawa sector of the electroweak standard 
model is the non-perturbative determination of the A-dependence of the upper and lower bounds 
of the Higgs boson mass |2p as well as its decay properties, where A denotes the cut-off of the 
theory. There are two main developments which warrant to reconsider these questions: first, with 
the advent of the LHC, we are to expect that properties of the standard model Higgs boson, such 
as the mass and the decay width, will be revealed experimentally. Second, there is, in contrast 
to the situation of earlier investigations of lattice Higgs-Yukawa models [||, ^ ^, ^, a consistent 
formulation of an Higgs-Yukawa model with an exact lattice chiral symmetry based on the 
Ginsparg -Wilson relation [||], which allows to establish a lattice version of chiral symmetry while 
lifting the unwanted fermion doublers at the same time. 

Before addressing the questions of the Higgs mass bounds and decay properties, we started 
with an investigation of the phase structure of the model in order to obtain first information about 
the region of the (bare) couplings in parameter space where eventual simulations of phenomeno- 
logical interest should be performed. 

In the present paper we basically summarize some of the most important results of our work 
on the model's phase structure, which we have studied analytically in the large A/y-limit for small as 
well as for large values of the Yukawa coupling constant [^], and numerically by means of HMC- 
simulations [[l0|]. Finally, we give a brief outlook towards some first and very preliminary results 
on the upper Higgs boson mass obtained at one selected cut-off A. 



2. The model and its numerical treatment 



The model, we consider here, is a four-dimensional, chirally invariant SU (2)l x SU (2)r Higgs- 
Yukawa model discretized on a finite lattice with L lattice sites per dimension. The model contains 
one four-component, real Higgs field and Nf fermion doublets represented by eight-component 
spinors \jf^'\ i = \, ...,Nf with the total action being decomposed into the Higgs action 5$, 
and the fermion action Sf. It should be stressed here that no gauge fields are included within this 
model. 

The fermion action Sp is based on the Neuberger overlap operator ^("^') [11] and can be written 

as 



Sf 



Nf 



(0 



2p 



('■) 



B„ 



(2.1) 

It describes the propagation of the fermion fields as well as their coupling to the Higgs field 
through the Yukawa coupling matrix B„ ,„ and the Yukawa coupling constant y^. Here the Higgs 
field <!>„ was rewritten as a quatemionic, 2x2 matrix 0„ = <I>°1 — i^irj, with f denoting the vector 
of Pauli matrices, acting on the SU (2) index of the fermionic doublets. 

Note that in absence of gauge fields the Neuberger Dirac operator can be trivially constructed 
in momentum space, since for that case its eigenvalues v^{p), e = ±1 for the allowed four- 
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component momenta p £ ^ are explicitly known. This will be exploited in the numerical con- 
struction of the overlap operator. 

The model then obeys an exact, but lattice modified, chiral symmetry according to 



(0 



1 



t 



-2/£0^ (2.2) 



recovering the chiral symmetry in the continuum limit 

The lattice Higgs action is given by the usual lattice notation 

5o = -K5v£<I>j, i^n+H+^n-il] + £<I>I<I>„ + A;v£ {^l^n -Nff 
n,n n n 



(2.3) 



with the only particularity that the fermion generation number Nf appears in the quartic coupling 
term which is a convenient convention for the large Nf analysis. However, this version of the lattice 
Higgs action is equivalent to the usual continuum notation jro|]. 

For the numerical treatment of the model we have implemented an Hybrid-Monte-Carlo (HMC) 
algorithm for even values of Nf with Nf/2 complex pseudo-fermionic fields (Oj according to the 
HMC -Hamiltonian 



Nfll 



7=1 



-1 



(2.4) 



where t, denotes the real momenta, conjugate to the Higgs field Since we focus here on checking 
the validity of our analytical investigation of the phase structure, which was determined in the large 
A'z-limit, the restriction to even Nf does no harm. For the further details of this HMC algorithm we 



refer the interested reader to Ref. [10|. 



The observables we will be using for exploring the phase structure are the magnetization m 
and the staggered magnetization s. 



m ■■ 
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(2.5) 



and the corresponding susceptibihties Xm = L''^ ■ [("J^) — ("^)^] and Xs = ■ [{^^) — i^)^] , where 
(...) denotes the average over the O-field configurations generated in the Monte-Carlo process. 

To locate the phase transition points, we decided to fit the data for the susceptibilities Xm> Xs 
as a function of according to the - partly phenomenologically motivated - ansatz 



Xn 



1 



7/2 



K 



(2.6) 



cnt ' 



where A™'*, 3 , and k^'j* are the fitting parameters for the magnetic susceptibility and staggered 
susceptibility, respectively, and v, 7 denote the critical exponents of the O^-theory. Here A™ 3 
(AI3) is actually meant to refer to two parameters, namely A™ (A|) for K^v < k"^^^ (Kn < ^cnt) 
Aj* (A3) in the other case, such that the resulting curve is not necessarily symmetric. The phase 



transition point is then given at the value of K^v 
susceptibility develops its maximum. 



K'crit i^N = K"crit) where the magnetic (staggered) 
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3. Large A//-limit for small Yukawa coupling parameters 

The phase structure of the considered Higgs-Yukawa model can be accessed in the large A^j- 
Umit by scaling the coupling constants and the Higgs field itself according to 

3^^ = ^^, ^N=T^, KN = kN, ^n = ^/N'f^n, (3-1) 

where the quantities yi^, Xn, K5v, and <i>„ are kept constant in the limit Nf — > oo allowing to factorize 
the fermion generation number Nf out of the effective action ^^//[O] = — A^/logdet(.^). 

One is thus left with the problem of finding the absolute minima of ^^//[O] in terms of the 
latter quantities. For sufficiently small values of the Yukawa and quartic coupling constants the 
kinetic term of the Higgs action becomes dominant allowing to restrict the search for the absolute 
minima of ^^^^^[O] to the ansatz 

^„=^- (^m + S-i-l)^""^ , OgR"^, |4>| = 1 (3.2) 

taking only a magnetization m and a staggered magnetization s into account. After some work, 
which was presented in detail in [^], one finally finds for the effective action 

Seffm = -NfZ log 



V+(P)| • \v+ip)\+J^{m^-f) ■ \v+ip)-2p\- |v+(p) -2p|) 



i 

while the Higgs action in this setting reads 



+ S<t,[^], (3.3) 



S.f =NfL'^-{ -SkN(m'^-f)+m^ + f + lN(m'^ + S* + 6m^s'^-2{m^ + f)) }. (3.4) 



The resulting phase structure in the large A/j-limit can then be determined by minimizing the 
effective action with respect to m and s. It is presented in Fig. |I]a for the selected value of the quartic 
coupling constant Xn = 0.1 and L = oo. Here we distinguish between the following four phases: 
(I) Symmetric (SYM): m = 0,S = (II) Ferromagnetic (FM): m^O,S = 

(III) Anti-ferromagnetic (AFM): m = 0,s^O (IV) Ferrimagnetic (FI): m / 0, 5" / 

In Fig. |T]b we compare this analytically obtained Nf = oo, L = oo phase structure with the re- 
sults of corresponding HMC-simulations performed on 8"*- and 6'^-lattices at Nf = 10. As expected 
we observe a good qualitative agreement between the numerical and analytical results. On a quan- 
titative level, however, the encountered deviations in Fig. |l|b need to be further addressed. These 
deviations can be ascribed to finite volume effects as well as finite Nf corrections. 

The finite size effects are illustrated in Fig. ^a, showing some phase transition points from the 
FM to the SYM phase as obtained from our numerical simulations on a 4^-lattice (open squares), 
and on an 8^-lattice (open circles) for the (very large) value of fermion generations Nf = 50, chosen 
to isolate the finite size effects from the l/Nf corrections. One clearly observes that the phase 
transition line is shifted towards smaller values of the hopping parameter when the lattice size is 
increased. The numerical results are compared to the Nf = oo phase transition lines obtained for 
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(a) (b) 

Figure 1: Phase diagrams with respect to the Yukawa coupUng constant and the hopping parameter K^v 
for the constant quartic coupHng Xn = 0.1. The black soUd Hne indicates a first order phase transition, while 
the remaining transitions are of second order [p^. (a) Analytically obtained phase diagrams for L — oo and 
Nf = °°. (b) Comparison with numerically obtained phase transition points for Nf = 10 and L'^ = 8^ (open 
squares) and = 6'^ (open circles). 




(a) (b) 

Figure 2: (a) Some selected phase transition points between the ferromagnetic and the symmetric phase, as 
obtained at Nf = 50 on a 4'*-lattice (open squares) and on an S'^-lattice (open circles), are compared to the 
L = 4 (dotted), L = 8 (dashed), and L = °o (solid) phase transition lines determined analytically in the large 
Nf -limit, (b) A^/^-dependence of k'^^^, k^^^ at = 2.0 as obtained on an 8^-lattice (open squares) and on a 
6^-lattice (open circles). The analytical, finite volume, large Nf predictions for the SYM-FM (SYM-AFM) 
phase transitions are represented by the dashed (dotted) lines. The dash-dotted lines are fits of the numerical 
data to linear functions as explained in the main text. In both plots = 0.1 was chosen. 



L = 4 (dotted line), L = 8 (dashed line), and L = oo (solid line). These analytically obtained lines 
perfectly describe the numerical results and one clearly observes the convergence of the numerical 
results to the analytically predicted L = oo line as the lattice size increases. 

The A^/-dependence of the numerically obtained critical hopping parameters and K^^-^^ is 
shown in Fig. ^ for = 2 . One clearly sees that for increasing Nf the numerical results converge 
very well to the analytical finite volume predictions, as expected. It is interesting to note that the 
leading term in the finite Nf corrections, i.e. the l/Nf contribution, seems to be the only relevant 
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correction here, even at the small value Nf = 2, as can be seen in Fig. ^ by fitting the deviations 
to the function fm,s{Nf) =A,„_s/Nf with Am,.v being the only free parameter Furthermore, one 
observes that the critical hopping parameter K^^^ is shifted towards larger values with decreasing 
Nf while K"^j.jj is shifted towards smaller values. 

For an investigation of the model at large values of the Yukawa coupling constant see Refs. [^, 

113. 



4. Outlook towards Higgs mass bounds 

In contrast to the previous discussion, where we considered the model mostly in the large Nf- 
limit, we now turn towards the physically interesting situation Nf = I. In order to investigate the 
model also at odd values of Nf we have implemented a PHMC-algorithm, which we will discuss 
in detail in an upcoming publication. 

The main goal here is to compute the cutoff A-dependence of the Higgs boson mass by fixing 
the top quark mass and the vacuum expectation value v to their phenomenologically known values, 
i.e. mtop = 175 GeV and v = 246 GeV. From this dependence one can eventually determine an 
upper bound of the Higgs boson mass. The v measured on the lattice has to be renormalized by the 
Goldstone renormalization factor Zq which can be obtained from the Goldstone-propagator G{p^) 
according to 

G-\p'') = ^ (4.1) 
Zg 

with denoting the squared lattice momentum. For the chosen setting {Kj^ = 0.240, = 0.711, = 
1 .0) we obtain Zg = 0.9662 it 0.0001 from the inverse Goldstone-propagator, plotted in Fig. ^a, and 
A = ( 1 684 ± 14) GeV. In Fig. lb we show one selected component of the fermion correlator {'^ft^ Yti ) 
yielding the top mass mtop = (170 ± 6) GeV in accordance with the phenomenological value. 




(a) (b) 

Figure 3: (a) Inverse Goldstone propagator G^'(p^) versus the squared lattice momentum fitted to a 
linear function, (b) Fermion time slice correlator {'^f,^^rf^) versus distance in time direction Af = \t2 — fi| 
fitted to a coi/i-function. 



In the presented setup we chose the relatively large value of the quartic coupling constant 
Aa? = 1 , aiming for an upper Higgs mass bound. In Fig. ^a we present the corresponding result for 
the Higgs correlator {f^n^t^) versus A?. We determine the Higgs mass by calculating the effective 
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mass at several values of At and finding its plateau value as shown in Fig. From this setup 
we find mu = (565 ±15) GeV. 

However, we remark that here we give only a first and very preliminary result towards our goal 
mentioned above. In particular, the value L ■ mtop = 1.62 is too small to determine the top quark 
mass reliably. Furthermore, the statistics (2500 configurations for the Higgs analysis) is still to low 
to obtain sufficiendy precise results for the physical quantities of interest. 
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Figure 4: (a) Higgs time slice correlator (O,, <I>,2 ) versus Af fitted to a coi/i-function. (b) Effective masses 
at Af fitted to plateau value nia. 
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